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We analyze and compare three diflterent strategies, all aimed at controlling and eventually halting 
decoherence. The first strategy hinges upon the quantum Zeno effect, the second makes use of 
, frequent unitary interruptions ("bang-bang" pulses and their generalization, quantum dynamical 

' decoupling), and the third of a strong, continuous coupling. Decoherence is shown to be suppressed 

f"**) ' only if the frequency A'^ of the measurements/pulses is large enough or if the coupling K is sufficiently 

, strong. Otherwise, if A or A' are large, but not extremely large, all these control procedures 

, ^ ■ accelerate decoherence. We investigate the problem in a general setting and then consider some 

' practical examples, relevant for quantum computation. 

PACS numbers: 03.67.Pp; 03.65.Xp, 03.65.Yz; 03.67.Lx 
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INTRODUCTION 



Interactions with the environment deteriorate the purity of quantum states. This general phenomenon, known as 
decoherence 0], is a serious obstacle against the preservation of quantum superpositions and entanglement over long 

■ periods of time. Decoherence entails non-unitary evolutions, with serious consequences, like a loss of information 
' and/or probability leakage towards the environment. 

, This issue is recently attracting much attention in view of interesting applications: for instance, the possibihty of 

■ controlling and eventually halting decoherence is a key problem in quantum computation 0] , where several computa- 
, tional states are simultaneously described by a single wave function and parallel information processing is carried out 

■ by unitary operations. In such a situation, efficient quantum algorithms need large scale computations, performed 
Qh' over (microscopically) long time spans Q. 

^Jj ■ A number of interesting schemes have been proposed during the last few years in order to counter the effects of 
Ch ] decoherence. Among these, there are quantum error correcting codes 0, schemes based on feedback or stochastic 
^ control ^ , the use of decoherence free subspaces and noiseless subsystems and mechanisms based on "bang-bang" 
qh (BB) pulses and their generalization, quantum dynamical decouphng 0,IEMEiEl- In this context, it was recently 
• • , proposed [T^ that the method of dynamical decoupling can be unified with the basic ideas underlying the quantum 

_ > ■ Zeno effect (QZE) [H Q (for a review, see llSj. |l6l). In particular, the decoherence free subspace is one of the 
, dynamically generated quantum Zeno subs p ac es within which the dynamics is not trivial (18J and whose subtle 

■ mathematical aspects are still debated [T^ll9ll20ll2ll l2^l23|. 

I It is worth stressing that the "bang-bang" scheme is a well-established "classical" control method, typically used 
in engineering problems and in connection with spin-echo techniques: see, for instance, Ref. p3 |. Its revival in 
quantum- information-related problems is only very recent. The key ingredient of BB and dynamical decoupling is to 
apply frequent (unitary) interruptions during the evolution of the system, in order to suppress the system-environment 
interaction. The similarity with the QZE was already noticed in It is however clear that the two procedures are 
physically equivalent, if one adheres to the commonly accepted interpretation of the QZE as a bona fide dynamical 
process, that can be completely explained in terms of unitary evolutions |2^. One should notice that this idea hinges 
upon a seminal remark by Wigner |2^, who introduced in 1963 the notion of "spectral decomposition," namely a 
dynamical process that associates a different wave packet with each eigenvalue of the observable to be measured. For 
example, the interesting proposal by Cook and the subsequent experiment with Rabi oscillations can be easily 
interpreted in fully dynamical terms when one observes that the "measurement" was realized as a dynamical process 
(optical pulse irradiation) [3> HE US ■ 

Once this physical equivalence is appreciated, the next logical step is a natural one: after having analyzed and 
understood the consequences of frequent unitary pulses, one studies the effect of a strong (unitary) continuous coupling. 
The relationship between these two procedures can be made mathematically precise (see Sec. ]^ and is of interest 
in itself: if an external field or "appar atus" is coupled to the system in such a way that the state of the system is 
"monitored" in some sense [sE IS l32L Is^ , a Zeno-like dynamics takes place in the strong coupling limit and once 
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again one can tailor decoherence-free subspaces This happens to be one of the most efficient and convenient 

control procedures, from a practical point of view. 

The aim of this article is to investigate these different physical procedures (Zeno, BB-dynamical decoupling and 
continuous coupling) and compare their effects. We will study the dynamics generated by very frequent interruptions 
(projective measurements or unitary "kicks," yielding dynamical decoupling), or by very strong coupling, and inves- 
tigate the possibility of designing decoherence-free subspaces. The method is general and can be applied to diverse 
situations of practical interest, such as atoms and ions in cavities, organic molecules, quantum dots and Josephson 
junctions (Mji35j 36, 3^1 . 

Our main objective is to endeavor to understand whether it is possible to control decoherence Is^ liolEll li ^ . 
Clearly, this requires a thorough understanding of the physical mechanisms that provoke decoherence and in general 
dissipative phenomena. One finds that very frequent kicks/measurements or very strong couplings can indeed control 
the evolution of the system and suppress decoherence. The physical mechanisms at the origin of this phenomenon are 
very close to the quantum Zeno effect. However, if the kicks/measurements are not extremely frequent or the coupling 
not extremely strong, both controls may accelerate decoherence. This extends the notion of "inverse" quantum Zeno 
effect (IZE) 43, '44] to a wider framework (not necessarily based on projection operators and non-unitary dynamics) 
and entails a deterioration of the performance of these schemes. We will analyze this effect in great detail and see 
that in order to avoid it, one must carefully design the control and study the timescales involved. Our analysis is of 
|!;en eral validity; however, for the sake of definiteness, we will study in particular the control of thermal decoherence 

This article is organized as follows. In Sec^Jwe briefly review the main features of the different control procedures. 
Our analysis is based on a master equation which is derived in Sec. liiil where the relevant timescales are emphasized 
and the general type of interaction specified. We then consider the case of thermal decoherence, discussing the Zeno 
control, the control via dynamical decoupling and the control by means of a strong continuous coupling in Sees. IIVI 
IVland lVll respectively. Some relevant examples are then considered in Sec. IVIll where we focus on the primary role 
of the form factors of the interaction in order to compare the different control procedures. Section IVlIII is devoted to 
conclusions and perspectives. 



II. CONTROL PROCEDURES: GENERALITIES 



Let the total system consist of a target system and a reservoir and its Hilbert space Htot — 'Hs (8i Tis be expressed 
as the tensor product of the system Hilbert space Tis and the reservoir Hilbert space TLb- The total Hamiltonian 

Htot^Ha+HsB = Hs®tB + '^s®HB+HsB (1) 

is the sum of the system Hamiltonian Hs (E" Is, the reservoir Hamiltonian I5 (g) Hb and their interaction Hsb, which 
is responsible for decoherence; the operators Is and 1b are the identity operators in the Hilbert spaces TLs and TYs, 
respectively, and the operators Hs and Hb act on Tis and T-Lb, respectively. 

The dynamics of the total system is conveniently expressed in terms of the Liouville operator (Liouvillian) >Ctot, 
defined by 

CtotP = -i[Htot,p] = ~i {HtotP ~ pHtot) , (2) 
where p is the density matrix. If the Hamiltonian is given by 1^, the Liouvillian is accordingly decomposed into 

^tot = Ca+ CsB ~ Cs + + CsBi (3) 

where the meaning of the symbols is obvious. We will not explicitly write the coupling constant A multiplying the 
interaction Liouvillian Csb- 

We focus on a proper subspace Ticomp C "Hs, in which quantum computation is to be performed. For this reason 
we will look in detail at the case 

'Hs = Hcomp ffi 'Worth- (4) 

In particular, when we will look at some concrete examples, in Sec. IVIll the computation subspace will be a qubit, 

^comp — ^ • 

Since, in general, the reservoir state is mixed, it is convenient to describe the time evolution in terms of density 
matrices. In the case of a quantum state manipulation, the initial state of the total system p{{)) is set to be a tensor 
product of the system initial state cr(0) and a reservoir (usually equilibrium) state pB 

p{Q)=cj{Q)®pB. (5) 
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The validity of this assumption, usually taken for granted, is discussed in Appendix ^ (see also 01)- The system 
state a{t) at time t is given by the partial trace of the state p{t) of the whole system with respect to the reservoir 
degrees of freedom: 

a{t) = iTBp{t). (6) 

When a{t) is not unitarily equivalent to ct(0) for a given class of initial states, decoherence is said to occur. The 
purpose of the control is to suppress such decoherence. Note that, for the control of decoherence, it is not necessary 
to look at all possible states: rather, it is sufficient to consider only those initial states which are relevant to the 
quantum state manipulation in question. 

A. Quantum Zeno control 

We first look at the Zeno control, by adapting the argument of Ref. (17|. The control is obtained by performing 
frequent measurements of the system. The measurement is described by a projection superoperator P acting on the 
density matrix 

p^Pp = ^(P„® 1b)p(F„® Is), (7) 

n 

where {P„} is a set of orthogonal projection operators acting on Tis- In the following, we restrict our analysis to a 
measuring apparatus that does not "select" the different outcomes (nonselective measurement) with a complete 
set of projection operators P„ — I5. The measurement is designed so that 

PHsB = ^ (Pn ® Is) HsB [Pn ® Is) = 0. (8) 

n 

In terms of the Liouvillian, this condition reads 

PCsbP = 0. (9) 

(We will see in the next subsection that a similar requirement is necessary for the BB control and for the control via 
a continuous coupling.) The Zeno control consists in performing repeated nonselective measurements at times t = kr 
[k = 0, 1, 2, . . .) (we include an initial "state preparation" &i t = 0). Between successive measurements, the system 
evolves via i?tot- The density matrix after A^ + 1 measurements, with an initial state p(0), is given by 

p{t) = p{Nt) = (Pe^'°'^pY p{0) . (10) 
We take the limit r ^ while keeping t = Nt constant and get 

p{t) =p[l + PCtotPr + O (t2)] * p(0) Pe^^'-^V(O) ■ (11) 

Equation © yields PCtotPp = -iP[Htot,Pp] = -iP[(Pi?tot), p] = -iP[H's Is + Is ^s, p] , whence 



PgP£totPt^(0) = Pe^'='V(0) 

= p(e-^-f^st0e-*^«*p(O)e*-^s*,g,e*^«*) , (12) 
where the controlled Hamiltonian -ffj'ot and Liouvillian are given by 

H[^, = PHtot = (Pn ® Is) i^tot (Pn ®tB)^H's®lB + ts®HB , (13) 

n 

C[^, = PCtotP = PCsP + CbP = C's + CbP ■ (14) 

with H'g = PHg = PnHsPn- Hcnce, as a result of infinitely frequent measurements, the system- reservoir coupling 
is eliminated and, thus, decoherence is halted. We notice the formation of invariant Zeno subspaces |l7j |: in the limit 
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of very frequent measurements, the evolution is given by (|13|l - (|14|) and transitions among different sectors of the 
Hilbert space vanish, yielding a superselection rule. The subspaces are defined by the superoperator P defining the 
measurement. The "decoherence-free" subspace is one of these Zeno subspaces. 
We will assume for simplicity that P commutes with the system Liouvillian 

PCs - CsP, (15) 

i.e. Hg = PHs ~ Hs, because our purpose is to control decoherence and we are not interested in a QZE over the 
system Hamiltonian Hs- The above assumption is equivalent to the following hypothesis on the Hamiltonian 

[Pn,Hs]^0, V?i. (16) 

In such a case 

= i^s + Cb)P . (17) 



B. Control via Quantum Dynamical Decoupling and "Bang-Bang" Pulses 

We now turn our attention to the so-called quantum dynamical decoupling 0, ^| , of which "bang-bang" pulses 
can be viewed as a particular case. The control of decoherence is achieved via a time-dependent system Hamiltonian 
He{t): 

H{t)=Htot+H,{t)(E)lB , (18) 
where Hc{t) is designed so that Uc{t) = Texp (^—i Hc{s)ds^ (T denotes time ordering) satisfies 

Ucit + t) = Uc{t) (19) 
dt{uUt) ® tB)HsB{Uc{t) ® Is) = 0. (20) 



In the interaction picture in which Hc{t) is unperturbed, the density matrix at time t = Nt, with initial state p(0), 
is given by p{t) = UtotiNT)piO)Ul^iNT) where 



UtotiNr) =T exp (-ij^ i?tot(s)rfsj = T exp ^-i^ i7tot(s)ds^ 



1 N 



(21) 



and i/tot(i) — {U^{t)^'lB)Htot{Uc(t)^'ll-B)- The second equality follows from the periodicity of i7tot(^)- A standard 
Magnus expansion of the time-ordered exponential jl^ leads to 



Texp ^-i Ht,tis)dsj = e-'^"' (22) 

where S'^^^ = i /J" Htot{s)ds and the term is of order (j = 1, 2, • • •). By assumption H20fl . one has 

i/(0) =H's®tB + ts®HB^ i/t'ot , (23) 

which is formally identical to ((T^ . where Hg = i dtUl{t)HsUc{t) = dxU^{xT)HsUc{xT) is independent of r 
because Uc{t) is r-periodic by H19|) and is always written as a function of t/r: Uc{t) = V{t/T). Therefore, in the limit 
T ^ while keeping t — Nt constant, one obtains 

C/tot(t) = [1 - ^Hi^T + O (r2)] ^ e-'^-t* = e"^^^* ^ e-'"^* . (24) 

In short, as a result of the infinitely fast control, the system-reservoir coupling is eliminated and, thus, decoherence 
is halted. As we shall see in a while, this is a consequence of the formation of invariant subspaces. 

As is well known, dynamical decoupling is a generalization of the evolution obtained by acting on the system with 
"bang-bang" pulses [Sj. In the latter, particular case, one applies during a time interval t two instantaneous unitary 
operators C/k and Ul and gets 

Hi, = Mot - ® i^tot [Pn ® Is) (25) 
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[see H13|l ]. where the projections P„ arise from the spectral decomposition 

[/k = ^ e-'^"P„, (A„ ^ Xm mod 27r, for n ^ m) . (26) 

n 

Notice that the map P is in this case the projection onto the commutant 

Z([/k) = {X| [X,C/k] =0}. (27) 

Equation l|25(l yields a convenient explicit expression of the effective Hamiltonian. As in the case discussed in the 
previous subsection, one observes the formation of invariant Zeno subspaces: transitions among different subspaces 
vanish in the r ^ limit, yielding a superselection rule. In this case, the subspaces are defined by (|25() - (|26|l and are 
nothing but the ergodic sectors of [/k- 

By assuming again, as in H15|) - (|16|l . that PHs = Hs and that PHsb = 0, as in JSJ-linj), we get that the controlled 
evolution for r ^ is given by 

Utot{t) = e-'"^"^* ^ e-^'"^* ® e-*^«* (28) 

or, in terms of Liouvillians, by 6^'°'* with C[^^ = PCtotP = {Cs + ^b)P, exactly as in l(T7|) . 

Moreover, in Ref. it was shown that one can obtain the same result H25|) by repeating a single "bang" , i.e. by 
using a single instantaneous unitary operator U^, without closing the group with ul- For simplicity, in the following 
we will always consider such a situation and will assume the commutation property IjlSI) . In such a case, the evolution 
is conveniently expressed in terms of the Liouvillian and density matrix 

pit) = [e^-e^'"'^] ^ Pp{0) e^'°'*Pp(0), r ^ , (29) 

where £k is the Liouvillian corresponding to the evolution H2t)|l and is given by H17|l . Note that the controlled 
Hamiltonians for bang bang pulses, (|25|l . and for the Zeno control, H13|l . coincide when the set of orthogonal projections 
{Tj) is chosen equal to the set (|26|l of eigenprojections of f/k, namely 

UP = 0, (PI) = 1. (30) 

Therefore, the two controls are equivalent in the ideal (limiting) case [T^ . However, throughout this article, by 
dynamical decoupling we will refer to a situation where the evolution is coherent (unitary), while by Zeno control to 
a situation where the evolution involves incoherent (nonunitary) processes, such as quantum measurements. 

The index "k" in the above expressions stands for "kicks." In the following, we shall use the expressions "bang-bang" 
pulses and "kicks" interchangeably. The latter is reminiscent of quantum chaos "4^. In fact, there is an interesting 
link between quantum chaotic dynamics, quantum diffusion processes and (inverse) quantum Zeno effect (501] . We will 
not elaborate on this issue in the present article. 

C. Control via a strong continuous coupling 

The formulation in the preceding sections hinges upon instantaneous processes, that can be unitary or nonunitary. 
However, as explained in the Introduction, the basic features of the QZE can be obtained by making use of a continuous 
coupling, when the external system takes a sort of steady "gaze" at the system of interest. The mathematical 
formulation of this idea is contained in a theorem |17| on the (large- iiT) dynamical evolution governed by a generic 
Hamiltonian of the type 

HK^Htot + KH,®tB. (31) 

which again need not describe a bona fide measurement process: can be viewed as an "additional" interaction 
Hamiltonian performing the "measurement" and A" is a coupling constant. 
Consider the time evolution operator 

UkU) = exp{-iHKt). (32) 

In the infinitely strong measurement (infinitely quick detector) limit K —>■ oo, the dominant contribution is 
eyip{—iKHct). One therefore considers the limiting evolution operator 



U{t)^ lim exp{iKHct)UK{t), 

K~^oo 



(33) 
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that can be shown to have the form 

U{t) - exp(-zff;„ti), (34) 

where 

ift'ot = Mot = Y.^Pn ® tB)Htot{Pn ® Ifi) , (35) 
n 

Pn being the eigenprojection of He belonging to the eigenvahie r]n 

Hc = ^ VnPn, iVn Vm, for fi ^ m) . (36) 

n 

By designing He so that PHsb = 0, the system-reservoir couphng is ehniinated and, thus, decoherence is halted. 
Equation 135|) . restricted to the system of interest, is formally identical to ()25|l and ()13f) . 
In conclusion, the limiting evolution operator is 



UK{t) - eM-^KHet)U{t) = exp (^-iKtJ^VnPn «) Is - ^Pi^tot^ j ■ (37) 

The above statements can be proved by making use of the adiabatic theorem 51]. It is worth noting that the evolution 
in the strong coupling limit is known to force the system to "cling" to the eigenstates of the interaction '5^ . In this 
sense, one expects that the dynamics be dominated by He for K large. The above theorem clarifies how the structure of 
He determines the features of the dynamics. Once again, like in the two previous sections, one observes the formation 
of invariant Zeno subspaces, that are in this case the eigenspaces of the interaction H35|l - I|36|) : the block-diagonal 
structure of (|37|l is explicit. The links between the quantum Zeno effect and the notion of "continuous cou pling " to 
an external apparatus or environment has often been proposed in the literature of the last 25 years |30l l3ll l33l |53j . 
The novelty here lies in the gradual formation of the Zeno subspaces as K becomes increasingly large. In such a case, 
they are nothing but the adiabatic subspaces. In terms of the Liouvillian: 

p{t) = e(-^'^=+^'°')*Pp(0) ^ e^'='*Fp(0), K ^ oo , (38) 
[see ij^ ] where the notation is obvious and 

CeP = 0, (PI) = 1 (39) 

[see (|30|) ]. The Liouvillian C[eit — P^totP corresponds to PHtot — H'tot ^-^id, under the assumption H15() and is 
again given by lfT7|l . 



D. Controlled evolution and Zeno subspaces 



The three different procedures described in this section yield, by different physical mechanisms, the formation of 
invariant Zeno subspaces. This is shown in Fig.^ If one of these invariant subspaces is the "computational" subspace 
'Wcomp introduced in Eq. Q), the possibility arises of inhibiting decoherence in this subspace. 

Of course, in the r, K^^ limit, decoherence can be completely halted, according to Eqs. (fT^ - (fnjl . (P^ - (P5|l and 
(|35|l - (|36|l . However, the objective of our study is to understand how the limit is attained and analyze the deviations 
from the ideal situation. This will be done by studying the transition rates 7„ between different subspaces and in 
particular their r and K dependence (see Fig.0. We shall see that in general this dependence can be complicated, 
leading to enhancement of decoherence in some cases and suppression in other cases. For this reason, the physical 
meaning of the expressions r, ^ in this section must be scrutinized with great care. 



III. MASTER EQUATION 

We consider the time evolution when the initial state is factorized as in © and the reservoir equilibrium state has 
inverse temperature f3 

PB = ^exp(-/3iJB) = i^exp(-/3S^)|/i)(Ai|, {CbPb = 0) (40) 
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FIG. 1: The Zeno subspaces are formed when the frequency of measurements or BB pulses or the strength K of the 
continuous coupling tend to oo. The shaded region represents the "computational" subspace Ticomp C Ti.s defined in Eq. 
The transition rates 7„ depend on r or K. 



where Z = tr^e = ^ the normalization constant. We will assume throughout our analysis that the 

characteristic timescales of quantum state manipulation in the space Ticomp [see (0J1 are much longer than any other 
timescales, so that the process is well described by the van Hove "A^f limit 54, 55.i5^l57|. where A is the coupling 
constant between system and reservoir [see the comment after Eq. Q]. For instance, if we take the timescale of 
quantum state manipulation to be of order A~^ {~ to a Rabi period in Ticomp), then the other energies involved are 
at most 0(A°). We will look at some concrete examples in Sec. I VIII 

Following Gardiner and ZoUer 's^, we now quickly derive the master equation and set up our notation. The starting 
point is the decomposition of the Liouville equation with the aid of the projection operators 

Vp = tTB{p}® PB = a® PB, Q = l-r (41) 

where tr^ stands for the partial trace over the reservoir degrees of freedom and /Os is the equilibrium reservoir state 
(HHl . Note that = V and ^ q Moreover, 

PCs = CsV, VCb - CbV = (42) 

and we assume that 

VCsbV = 0, (43) 

which can always be satisfied by redefining the system Liouville operator Csp — > Csp + trB{CsBp} ^ Pb and the 
interaction Liouville operator LsbP — ^ ^sbP — ^t^bI^sbp} Pb- 
The evolution in the interaction picture reads 

Piit) = CsB{t)pi{t), CsB{t) = e-^«*£sBe^°* (44) 

and by applying the projection H41|l together with Eq. (|43() one gets 

Vpi{t) = VCsB{t)Qpi{t), Qpi{t) - QCsB{t)Vpi{t) + QCsB{t)Qpi{t). (45) 

By formally integrating the second equation and plugging the result into the first one, one obtains to order A^ 

Vpi(t)= [ dsVCsB{t)QCsB{.s)Vpi{s), (46) 
Jo 

where the initial condition jSJ, yielding Qpi(O) — 0, was used. By using the definitions (|41l) and the conditions 
(|23-(|1SI, Eq. gni) yields 

ai{t) = / dsJCi{t,s)ai{s), (47) 
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where 



/Ci(i, s) (T tvB {CsB{t)CsB{s) CF (g) pb} ■ 



(48) 



By making use of the first Markov approximation ai{s) <yi{t) [Sg, which is motivated by the fact that the bath 
correlation kernel /Ci(t, s) is different from zero only for s t — such that ai{t — Tc) ~ cTi(t), one gets 



&i{t) = C{t)ai{t), C{t)= dslCi{t,s) . (49) 

Jo 

If the time t in Eqs. (|49|l is much larger than the bath correlation time, t 3> Tc, one can safely replace the upper limit 
of integration with oo, getting a Markovian equation with the time independent Liouville operator C = C{oo). 
We emphasize that this procedure can be rigorously justified in the (weak coupling) limit [s^ 



C = lim 



(50) 



which physically corresponds to a time coarse-graining ansatz |59|,|60|. From 148|l and (|50|l one gets (by suppressing, 
for simplicity, the subscript I for the operators in the interaction picture) 



La ^ lim tr^ < e 



t/A2 



ds CsbCsb{s) 



trs S X! '^'^ j CsbCsb{s 



,a®PB 



(51) 



where are the eigenprojections of the Liouvillian £5, 



(52) 



and in the limit the off-diagonal terms, e*^'^ Q^[- ■ -IQu)' , vanish due to the Riemann-Lesbegue lemma. Notice 

that the superoperators Q^^ can be expressed in terms of the eigenprojections of the Hamiltonian as 



QmPQn-, Hs — EnQn 



(53) 



From a physical point of view, the result (|51|l hinges upon a second-order perturbation expansion of the Liouvillian 
(O in the interaction picture 



^ — Cot^Ctati 



= Texp( / ds £sb(s) J ~ 1 + / ds £sb{s) + ds dsi £sb{s)Csb{si)- 
\Jo / Jo Jo Jo 



(54) 



Indeed, the first-order term vanishes after the projection due to H43|l . while the projected second-order term reads 



trs ds J dsi Csb{s)Csb{si)cf ® pB^ = j 

-1° ° 



= ds tiB { e 



ds trs < e 





Ct a 



dsi CsB^SBisi) 



J dsi CsbCsb{si) e^^V^psj 
Pb| ~ i trs I ^ J ds CsbCsb{s) 



(55) 



In the second equality we considered times t much larger than the bath correlation time Tc, so that the integration range 
can be extended from (— s,0) to (— oo,0), while in the third equality we neglected the rapidly oscillating (compared 
with those responsible for decoherence) off-diagonal terms. By combining (|55|l and (|54() we finally get 



(n{t) =trB{e-^«*e^'°'*ai(O)0pB} ~exp(/:t)ai(0) 



(56) 



which is nothing but (|49|l . when one substitutes C{t) C{oo) = C. 

Some of these ideas and techniques, at different levels of rigor, have been investigated and applied in the literature 



of the last four decades I5J, l55|, l5 
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A. The general case 
Assume now that the interaction Hamiltonian Hsb can be written as [s^ 

HsB = ® ^™ + ® ^"0 ' 

m 

where the Xm are the eigenoperators of the system LiouviUian, satisfying 

CsXm = iuJmXrn (t^m ^ W„, for m ^ Tl) 

and Am are destruction operators of the bath 

Am = A{gra) = / (f' k g*m{k) a{k) , 



(57) 



(58) 



(59) 

expressed in terms of bosonic operators a(fe), with form factors grn{k)- We are specifying our analysis to three 
dimensions (although it is valid in any dimensions). Incidentally, the form of the Hamiltonian (|57|l is of very general 
validity (and is not limited, as one might naively think, to dipole-like approximations): the only assumption made is 
that the coupling with the bath be linear, i.e. one is not considering terms of the type a^, a^^, etc., which would only 
be relevant for squeezed reservoirs. In practice, one determines the operators (|58f) . then finds the bath operators in 
order to write the interaction in the form (|57|l . and neglects nonlinear terms. 
In the interaction representation we get 



HsB{t) = e-^«'HsB =Y,{Xrn® Al{t) + X ® A„(t)) , 

m 

where 

Am{t) = j d?k g:^{k) e-'^("'=-"-)*a(fc) . 
If the bath is in the thermal state (|40|l we obtain 

{Am{t)Al{s)) = j d^k |.9„(fc)p {N{uk) + 1) ".)(*-^) 
{Al{t)Am{s)) - / d'k |5,„(fc)|2Ar(^,)eH-.-™)(*-^) 



(60) 



(61) 



(62) 



and {Am(t)Am{,s)) = (AL(t)At„(.s)) = 0, with N{uj) = l/(e^^ - 1). 
From H51(l we get 



Ca = ds ti-B lQu^SB^SBis)QuCr Pb\ 



and by using the property 



Hi, 



which easily follows from the definition l|52(l . we get 

Ca^-J2 J ds ti B {[{q H s b) ,[{Q-u;Hs B{s)y a (S)pb]]} 
By using and (|57|) one obtains 



(63) 



(64) 



(65) 



^HsB = H^sB = X-m ^ Al^ + 



(66) 
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whence 

La — — / ds trs I 



(m) 
SB ' 



(-m) 
SB 



X 



l^® A^^{s),a® pB } .(67) 



In the second equahty we neglected terms containing two annihilation or creation operators, which identically vanish 
by performing the trace over the thermal state ps- Equation Ht)7|) can be put in the form |i58j 



Ca = ~i ^ [5ynX}^Xm + e„iX„iXl^, cr] 

m 



and 



^2 m ^^m 



dt{A^{Q)Al{t)) 



dt{Al{0)A„,{t)) 



(68) 



(69) 



The first line in H68|l is just the renormalization of the free Liouvillian Cs by Lamb and Stark shift terms and will be 
neglected in the following. The dissipative part is given by the second line, which appears in the Lindblad form, so 
that tr£cr = 0. 

In Eq. (|58|l we will identify cij_,„ = -^oJm, and will assume that X^m = and g,n — .9-m, which is equivalent to 
the hypothesis that the interaction Hamiltonian be the product of selfadjoint operators acting on the system and the 
bath, namely Hsb — Si -^s "* ® -^s^ with Hg^ and H'j^^ selfadjoint. Notice, therefore, that we are not making any 
rotating- wave approximation, and the interaction Hamiltonian Hs b H57|) contains both rotating and counter- rotating 
terms. The dissipative part of (|68|l can now be rewritten as 



(70) 



£cr = 70 ( XqoXo - - {XqXo, o-} 



^ 7m (x^'jXI - i {xt„X™, a}) + 7-™ {xlaX„, - \ {X,^Xl,a} 



where 7,„ = K„i + G_,„. We introduce the bare spectral density functions (form factors) 

'«m(t^) = j d^k\gm{k)\^5{uk-^^)i K„i(w) = 0, for cj < 
and the thermal spectral density functions. 



K^(w) = K™(tj) {N{u) + 1) + nr,^{~uo)N{-UJ) 



1 



1 - e-P'^ 

which extend along the whole real axis due to the counter-rotating terms and satisfy the KMS symmetry [5£ 

'«m(-^) = X , -, «m(^) = — TTT ^ '^ml^) = exp(-/3w) ^ (w). 



We explicitly get 



and 



N{uj) + l ' Ni-co) 

ra ) ~^ ^ } 1 



Gn 



7„ = 2Re / dt (^{A^Alit)) + (^Al^A_^{t)^^ = ^t^k^M 



(71) 



(72) 



(73) 



(74) 



(75) 
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FIG. 2: The form factors at zero temperature, K{Lo)/{g'^W) vs cu/W. Full line: exponential form factor H76|l : dashed line: 
polynomial form factor 1771 . 

It is useful to look at some concrete examples and scrutinize the modification of the form factor (|71|l . due to the 
presence of the thermal bath. Let us focus, for the sake of clarity, on two particular Ohmic cases: an exponential 
form factor 



4f H=g2^exp(~^/A)^H 



and a polynomial form factor 



UJ 



(76) 



(77) 



In the latter case, we focus on n = 2, which is typical of quantum dots [23 (the case n = 4 is also of interest, being 
the nonrelativistic form factor of the 2P-1S transition of the hydrogen atom IgiI, ^6^). In the above formulas, g is 
a coupling constant, A a cutoff and 9 the unit step function. In order to properly compare these two cases, we will 
require that the bandwidth be the same: 



W = 



cLj Km^ (Ul) diV Km' (w) 



(78) 



where the square root of the denominator 



duj Krn{uj) 



(79) 



is the so-called Zeno time, characterizing the convexity of the survival probability at the origin 1161 l62l l63| . Notice 
that a finite natural cutoff A ~ 8.498 • 10^*rad/s and a finite Zeno time tz — 3.593 • lO^^^s can also be computed for 
the hydrogen atom in vacuum [polynomial form factor H77|l with n = 4], as well as for atomic and molecular systems 
whose electronic wave functions are known. The condition 178|l when n = 2 yields the ratio Apoi/Acxp = 1.275 between 
the cutoffs for the polynomial and exponential form factors, and W = 1.99Aoxp. The two form factors are displayed 
in Fig.Hfor g = 1. 

The thermal form factors H72|l are displayed in Fig. |3 for two different temperatures. Three features are apparent: 
the form factor is an increasing function of the temperature /?~^. Its value at w = is k(^(0) = k^,j(0+)//9 — g'^/P, 
where the prime denotes derivative. Moreover, its derivative reads k^'(0*) = k^(0+)/2 ± kJ^(0+)/(2/?), whence it is 
continuous, ^^(0) = 5^/2, in the polynomial case (because k'^{0~^) = 0), and discontinuous, k(^'(0^) = g'^/2^g'^/ (/3A), 
in the exponential case; this is more apparent at higher temperatures. Finally, the support of the thermal form factors 
is no longer lower bounded, due to the effect of the counter-rotating terms. 
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FIG. 3; The thermal form factors ^^(aj) vs lj. FuU lines: exponential form factors 17211 . 17611 : dashed lines: polynomial 
form factors 1721 1. 1771 . The form factor is larger at higher temperature f3~^ . Note the discontinuity of the derivative in the 
exponential case at oi = (more apparent at high temperature). We set PW — 50 (low temperature) and l3W — 2 (very high 
temperature). 



B. Two- level system 



A particular case of the above is the qubit Hamihonian 



Ho = -a,, HsB ^cfz® (^(<?o) + A^{9o)) +<yx® (^(ffi) + ^^(ffi)) • 



This is of the form I57|) , when one identifies 



(80) 



Xo = CTz, X±i = (Tzp = , tJ±i = ±i2, Wo = 0, 



(81) 



hence 



Cp = 70 {(J^pa^ ~ P)+ 7+1 ( cr-PO-+ - ^ {a-+cr_, p} j + 7_i Ur+pd- - ^ {cr-cr+, p} ) , 



(82) 



with 



where we used (|7^ . 



^o(0+) 
ov"; =^ ^" — ^ — ; 



70 = 2t:k'q{Q) = 27r 



7±i = 2TiK'l{±Vt) 



IV. QUANTUM ZENO CONTROL 

Let us look at the quantum Zeno dynamics with a finite interval t = t/N between measurements, 



(83) 



p{t) = 



pe^tatrp 



P(0) 



(84) 



where £tot and P are given by 131 and 0, respectively. We will look at the subtle effects on the decay rate arising 
from the presence of the short-time quadratic (Zeno) region. Therefore the standard method (5^ is not applicable to 
the present situation and the limit must be evaluated by a different technique. We only sketch the main steps in the 
derivation and give more details in Appendix IBl Second order perturbation in £5^ and the conditions ljHJl-l|njl yield 



Pe 



Aotrp^ Pg£or-^gXp I / ds CsBis)] P ^e'''>''P 



Jo Jo Jo 



P . (85) 
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In terms of the operator 5z (t) , defined as the solution of the operator equation 



dse-^«^gz(T)e^«^ =P / ds ds, CsBis)CsB{si)P = ds e'^" 
Jo Jo Jo 



Cos 



dsi PCsB^SBi-Sl)P 



one obtains 



t/r 



Fe^o^Texp / ds e-''<>''gz{T)e''<>' P 



.t/r 



= Pexp{[/:o+az(T)]0 



(86) 



(87) 



Under the assumption that the bath state can well be approximated by an equilibrium state at time the final 
reduced state cr(t) is shown to satisfy the equation 



CT{t)^[Cs+Cz{T)]cj{t) 



with 



Cz{t)(j = tr_B {Qz{T)a (g) ps} 
Note that Cz{t) is the solution of the operator equation 



(89) 



dt : 



-Cst, 



Cz{T)e^''= / dtPL{t)P = / dt dsPICi{t,s)P. 



(90) 



'0 "'0 Jo Jo 

where JCi is defined in H48|l. The dissipative part of 1)89(1 is found to have the explicit form [analogous to Eq. 1)70(1 ] 



+ E ^™(^)^ (xmP'yXl - i [xlXr.,Pa]) + E 7^,„(t)P (xlPaX^ - i [x„,Xl,Pa}\ 

m>l ^ ^ m>l ^ ^ 



where the controlled decay rates read 

2 /-r /.O 

T 



7^(t) = ^Re^ dt^^rfs ((A™AL(s)) + (aL™A_„,(s))) = du; kUlu) sinc^ j^- 



(91) 



(92) 



with sinc(a;) = sin(x)/a;. This yields Zeno and inverse Zeno effects as r is changed, as we will see in Sec. lVIll The key 
issue, once again, is to understand how small r should be in order to get suppression (control) of decoherence (QZE), 
rather than its enhancement (IZE). 



V. CONTROL VIA DYNAMICAL DECOUPLING 

We can now investigate the nonideal bang-bang control of decoherence. From Eq. 1(29(1 . describing a BB control 
with a single kick |12ij. 

p(i)= [e^^e^-'^]^p(0) , (93) 

where £tot is again given by (O. As in the Zeno control, we consider here the case where t is finite, so that the 
effects on the decay rate arising from the presence of a short-time quadratic (Zeno) region play a fundamental role. 
Once again, we only sketch the main steps in the derivation and give more details in Appendix [CI Second order 
perturbation in Csb yields 

g£,g£totr ^ g£,g£or^exp ( /" ds Csb{s)) ^ c^'^e^"^ 1+f ds CsBis) + f dsf dsi CsBis)CsBisi) .(94) 

V^O / I Jo Jq Jo J 

In terms of the operators .^^(t) and G]^{t), defined as solutions of the operator equations 

/ dse-^-''T],{T)e^^' = f dsCsBis) (95) 
Jo Jo 



ds 



Sk(r)e 



C-rS 



ds / dsi 
Jo 



CsB{s)CsBisi) - e-^-'J^i,{T)e 



^k(T)e 



(96) 
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with 
one has 



Ct = h Co, 

T 



1 JV 



exp 



exp 



+ /:o+J^k(r)+e;k(r) 



With the aid of ()98fl . the final reduced state o'(i) satisfies the equation 



a{t) = 



a{t) 



with 



£k(T)cr = tl-B {C/k(T)CT (g) Pb} ■ 

The dissipative part of (|100|l has the exphcit form 

£k(T)a = ^1{t) fxo(r)aXo(r) - \ {Xo(r)Xo(r), a}) + ^ 7,^ (r) ('x„(r)aXt„(r) - ^ {Xt,(r)X„(r), a} 



(97) 



(98) 



(99) 



(100) 



E f^,UT)^^»(r) - i {X„,(r)Xt„(r)a}) , 



where, in analogy with Eq. H58|l . the Xm(r) are the eigenoperators of the Liouvillian C-^/t + £5, satisfying 

^ + £5 j X„i(t) = itj„(T)X„(r) {ujrn ^ ^n, for m ^ n) 



(101) 



(102) 



and the controlled decay rates read 
7,'j,(T)=2Re 

with 



r ((i,„(0)itjt)^ + (iL„.(0)i_„(<))) = 27r^^ (c.„(r)) = 2^«f„ + 0(1)) , (103) 



A^{t) = / d3fcg;^(fc)e-^(-^-"(-))*a(fc). 



(104) 



Notice that the mechanism of decoherence suppression (I103|l is not fully determined by £tot and P, in contrast to 
the Zeno case, and depends also on the details of the Liovillian £k through Li;„i(r). This is best clarified by explicitly 
looking at a particular case: let us consider the two level system (|Sn|) with 50 = (spin-flip decoherence). We include 
an additional third level — that performs the control — and add to (I80|) the Hamiltonian (acting on Tis ® span{|Af)}) 



Hm = -^\M){MI 



(105) 



so that \M) is degenerate with ||). The control consists of a sequence of 27r pulses between ||) and \M), given 
by 



where 



C/k = exp [-i^{\i){M\ + \M){i\)] = P^-P. 



Ft = I T)(T I, P-i = Pi+Pm = I i)a I + |Af)(M|, 



(106) 



(107) 



are the eigenprojections of f/k (belonging respectively to 6^*"^^ = 1 and e^^^^^ = — 1) which define two Zeno subspaces. 
In the r ^ limit any decoherence between these two subspaces is suppressed. In fact, the total decay rate of the 
upper level has been explicitly computed 0, l64j and reads 



71- (r) = lim t j duj k'^ iuj) sine 



t tan^ 



(108) 
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As a matter of fact, the function multiplying the thermal form factor inside the integral can be shown to have the 
interesting limit 



lim t sine 

t — ^oo 



tan^ - 



V 2 



1 



27r 



The above limit is taken by keeping r fixed — finite and nonvanishing — and t 
plugging H1U9|) into (|108|l one gets 



9 °° 1 



2tt 



[n + —ij + i/2)]+K,'^ in 



27r 



U + 1/2) 



Nt, with N integer and even 



27r 



(j + 1/2) 



(109) 



By 



(110) 



which is a sum of suitably weighted terms of the form l|103() . This yields again control of decoherence as r is varied, as 
we will see in Sec. lVIll The key issue, once again, is to understand how small r should be in order to get suppression of 
decoherence (control), rather than its enhancement. Equation (|110l) yields also a significant computational advantage, 
when compared to H108|) : for well-behaved form factors (without resonances) the first few terms already provide a 
good estimate of the controlled lifetime. 



VI. CONTROL VIA A STRONG CONTINUOUS COUPLING 



We can now analyze the last case, that of control by means of a strong continuous coupling. Since the control of 
decoherence is achieved by adding a control Hamiltonian KHc acting on the Hilbert space Tis, we begin with the 
study of the spectral properties of the new "system" Hamiltonian Hs{K) = Hs + KHc- By writing the spectral 
resolutions of Hs and H^ 



ifs-^£^nO„, Hc = Y,Vn.Pm, (HI) 

n m 

with J2n Qn = J2m ^™ ~ ^' ^^'^ using the property H16|) we see that Pmn = PmQn is a (finer) orthogonal resolution 
of the identity, i.e. „ Pmn = 1, with PmnPm'n' — ^m,m'5n,n'Pmn- Notc that some Pmn can vanish. In particular 
H'g can be explicitly diagonalized 

Hs{K) = ^ [KT^m + S„) Pmn. Pmn = PmQn ■ (112) 



Equations Hlll|) and ()112f) directly translate in terms of Liouvillian as 

£5^-1^ UnQn, £c ^ ^mPm (113) 

n ni 

and 

Cs{K) = Cs + KCc = -lY,^rnn{K)Pmn, 
UJmn{K) = K^m+i^n, Pmn = PmQn ■ (114) 

The condition ||SJ) for a complete control of decoherence, PHsb = 0, leads to 

= PHsB = J2 PmHssPm = PoHgB - PoQnHsB - PonHsB, (115) 

m n n 

whence 



PonHsB - 0, Vn . (116) 

Therefore, by following exactly the same steps of Sec. lIII with Hs{K) defined by 1112|) in place of 7?s, one obtains 
that the dissipative part of the Liouvillian Lk governing the slow evolution of the reduced density matrix a is given 
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by 



m>l,n ^ ^ 

+ l~-^AK)(xl^aX„,n-UXmnXl,^,c7}) , (117) 



771 > 1 , n 



where 



^m77 = P-mXn , (118) 

with X„ given by lt57|) - (|SS|l . and 

-l^n{K) = 2Re^ ((i,„„(0)ij„„(0) + (iL„,„(0)i_™,„(0)) = 2^Kf (w™„(i^)) , (119) 

with 

A^mn{t) = I d3fc5,:(fc)e-'("--"'""(^))*a(fc) . (120) 



All terms with to = identically vanish due to pi6f) . In the K ^ oo limit, because the thermal form factor k(^(w) 
vanishes as lo oo (cf. Fig. 2), one has 

-1ran{K) = 2^K^ {Kilm + UJn) ~ 17^^ (ifl^™) , for OO . (121) 

Hence, in the K —>■ +00 limit, the dissipative part disappears, Ck 0, or decoherence is suppressed, as expected. 

It is interesting to observe that, when the condition (|16|l is not satisfied, the control via a strong continuous coupling 
needs an additional argument. In such a case, the control Hamiltonian He and the system Hamiltonian Hs cannot be 
simultaneously diagonalized, but (for a finite-dimensional Ti-s), as a result of the analyticity of the eigenvalues and the 
corresponding eigenprojections of the Hermitian operator Hs{K) / K = Hs/ K + He with respect to the perturbation 
parameter 1/K 65], the eigenvalues u)mn(,K) of the new system Liouvillian Cs{K) = KCc + Cs and the corresponding 
eigenprojections Pmn{K) satisfy 

UJmn{K) = K^m + ^^^i + O ( ^\ , (122) 



K 



P (K) - P^°'> + — + o( — ^ 



(123) 



where fimi and Pml (j = 0, 1) do not depend on K. As in H116|l . one gets that P^Hsb = 0, but this does not imply 
that Pon{K)HsB = 0. As a result, there appear dissipative terms which tend to via a different mechanism from the 
one outlined above. This aspect will be discussed elsewhere, together with similar phenomena that occur also for the 
other two control mechanisms (BB and Zeno). 

In general, as in the BB control but in contrast to the Zeno case, the mechanism of decoherence suppression (|121f) 
is not fully determined by Hs and depends on the details of the Hamiltonians Hs and He- Once again, this can be 
clarified by looking at a specific example: consider the two level system (|8U|I with go = (spin flip decoherence). We 
add to llSUj) the Hamiltonian (acting on Hs ® span{|M)}) 

Hm - ~^\M){M\+KH„ 
H, = |i)(M| + |M)(i| = P+ - P_ , (124) 

where 

P, = M±M±(M,|±)(±|. (125) 

The third state \M) is now "continuously" coupled to state ||), K £ M. being the strength of the coupling. As K is 
increased, state |Af) performs a better "continuous observation" of yielding the Zeno subspaces |l6| . In terms of 
its eigenprojections, H^ reads [see 

i7c = 77TA+'?-^-+'?+^+> (126) 
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with P| = |t)(TI Eind r^i = 0,rj± ~ ±1. In the Zeno hmit (_fC — > oo) the subspaces 7i|. 7i+ and H- decouple due to 
wildly oscillating phases 0{K). We get 

PHsB = P]HsbP^ + P-HsbP- + P+HsbP+ = 0. (127) 

Therefore in the limit K oo, j±i — and decoherence is halted. 
We can diagonalize the new system Hamiltonian 

H's = - §|M)(M| + KH^ ^ ^P, + ^-^ + i^) P+ + (-^ - P_ . (128) 

The new system operators (I81|) become 

X± = P±a,P| = -i=|±)(t|, Xo = \-){+\, C'sX±=2{nTK)X± , C'sXo^2tKXo, (129) 

HsB = {x+ +X^+Xl+ Xl) ® {A{g) + A^ig)) , (130) 

hence 

£kP = 1+{K) (^X+pXl - i +7-(i^) (^X-pXl - I {xlx^,p} 

+ ^+{K) (^xlpX+ - i [x+xIp]^ +1-{K) (^xipX^ - 1 , (131) 

where 

-/±{K) = 27rn^{nTK), j±{K) ^ 2TTK'^{-n ± K) . (132) 
For example, the decay rate out of state |t) reads 

VII. THE ROLE OF THE FORM FACTORS 

We can now test the general scheme described in the previous sections by looking in detail at some particular 
cases. We will consider the two-level situation and compare the three control methods both with exponential H7fci|) 
and polynomial form factors (|77|l . We will concentrate on the transition between a regime in which decoherence is 
partially suppressed ("controlled") and a regime in which it is enhanced. We shall work in a high-temperature regime, 
which is rather critical from an experimental point of view, because of temperature-induced transitions in two-level 
systems. We shall set ft = O.OIW and (3 = 50W^^, so that temperature^/3^^ = 2i}. 

A. Quantum Zeno control 

We first consider the Zeno control by projective measurements. Dissipation and decoherence are characterized by 
the decay rate (|^ : 

72(t) =T /°° d^Ac^(a;) sinc^ f^^i^r) ^ ^, r ^ 0, (134) 



where tz, 



dcj K^{uj) = / doj k((jj) coth I -— j ~ / duj k{uj) + 2 / duj k{uj) exp{— f3uj) , (135) 
; Jo V 2 / Jo Jo 
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FIG. 4: Projective measurements: ')'^{t)/^ vs Wt. Full line: exponential form factor 17611 : dashed line: polynomial form factor 
l|77|l with n = 2. 



is the thermal Zeno time. (We dropped the suffix m for simpHcity.) Observe that, by making use of the Umit 



r — yoo 

one gets 



hm T sinc^ Hr = 27r(5(w) (136) 



7^(t)^7, r^oo, (137) 

where 

7 = 27rK''(r2) (138) 

is the natural decay rate (|75|) . The ratio 7^(t)/7 is the key quantity: decoherence is suppressed if 7^(t) < 7, and 
it is enhanced otherwise. This ratio is shown in Fig. 2]as a function of t [in units H^-the bandwidth defined in Eq. 
(|75)l] . The transition between these two regimes takes place at r = r* , where r* is defined by the equation ^ij 

7Z(r*) =7^(00) =7. (139) 

If T* belongs to the linear region (|134|1 (which is our case and is true for sufficiently small energy Q. of the initial 

state), one gets 

r* ^ ^rl = 2n ^y\ (140) 

The short time region is displayed for clarity in Fig. |S1 

It is useful to spend a few words on the physical meaning of the expressions t — > 0, — > 00 in the above (and 
following) formulas. Times and temperatures are to be compared with the bandwidth W (or frequency cutoff A). 
Times (temperatures) are "small" when t ^ <C W). For example, it is worth emphasizing that the relevant 

timescale is r* , when one considers short-time expansions in a Zeno context |44l l63j : the expansion H134|) is valid for 
r < (and not r < rz, as it is sometimes erroneously assumed). 

B. "Bang bang" Control 

We now discuss BB. The decay rate is given by Eq. (|110|l : 

9 °° 1 

7V) = ^E7-^K(^ + >- + i))+-^^("-7(2^' + i)) 

j=Q [J + 2) 

j=o \J + 2) J=0 [J + 2) 
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FIG. 5: Projective measurements: 'y^{r)/'y vs Wr, for small r. Full line: exponential form factor 17611 : dashed line: polynomial 
form factor 1771 1 with n = 2. r* (indicated) is defined by the equation 7^ (r*)/7 = 1. Decoherence is suppressed when 7^ (r) < 7; 
it is enhanced otherwise. 



where we made use of (|73|l in the first expansion and assumed that (3 is not too small (as compared to t) in the second 
one. In the exponential case (|76|l one gets 

^{E) (l{2j + 1)) = g2^(2j + l)e-^(2j+i) ^ ^(E) + l)e-^^^ , (142) 

whence 

7V)--^^(^)(-), r^O, (143) 
while in the polynomial case H77|l one gets 



whence 



where C,{x) = ^^'^ Riemann zeta function. 

On the other hand, in both cases, 



7^ 



(144) 



4 °° 1 



j=o U + 2; 
where we summed the series 

°° 1 °° 1 TT^ 

V 9=4y ^ = 3C(2) = — . (147) 

The ratio 7'^(t)/7 is shown in Fig. as a function of r. Once again, the transition between the two regimes takes 
place at T = T* where r* is defined by the equation 

-f^{r*)=-/^{oo) = j. (148) 

If T* is in the asymptotic region (|143|l one gets in the exponential case (|76f) . 

(^) - = ^«^(^^), (149) 
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FIG. 6: BB kicks: 'y^{T)/'y vs Wt. Full line: exponential form factor 1176^ : dashed line: polynomial form factor l|77^ with n — 2. 
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FIG. 7: BB kicks: 7'^(t)/7 vs VKr for small r. Full line: exponential form factor I76II : dashed line: polynomial form factor I77II 
with n = 2. Decoherence is suppressed when 'y^ir) < 7; it is enhanced otherwise. 



which yields 



A V 8 g^A 



'T .g2A 



(150) 



where W is Lambert's VF-function j^ali that is the inverse of the function f{W) — We^ , and we have taken its — 1 
branch. 

On the other hand, for the polynomial case l(77|l one gets from (|145fl 



8(1 



^)C(2n+l)^ 4(1 



lC(2n+l) 



(151) 



and 



A V8(l - 2-2"-i)C(2n + 1)52A 



A V4(l - 2-2"-i)(;(27i+ 1) g'^A 



(152) 



The short-time region is shown in Fig. [7| It is useful to observe that the results H149|) - (|15U|I and H151|) - (|152|l bear an 
important dependence of r* on the "tail" of the form factor. This is to be sharply contrasted with the projective 
measurement situation (|140|) , that yields a dependence of the transition time r* on the "global" features of the form 
factor. This difference is apparent if one compares Figs. |31 and [7| and shows that the latter method offers important 
advantages if one aims at inhibiting decoherence, because of the larger (and easier to attain) value of r*. 
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FIG. 8: Continuous coupling: ^'^{K)/^ vs 2nW/K. Full line: exponential form factor 1)76^ : dashed line: polynomial form factor 
iTfTl with n = 2. 



C. Control by continuous coupling 



Finally, we can look at continuous coupling. The timescale for decoherence is (|133|l : 

/oc 
duj Kl^iuj) [6{uj - n - K) + 6{uj - n + K)] 

= 7r[K'^{n + K) + K^\n-K)] ^TTK{K){l + e-'^^) r~.TTK{K), K^oo. (153) 

On the other hand, 

-/''{K) ^7, K ^0. (154) 

Notice that the role of K in Eq. 11531) and the role of 1/r in Eqs. (|143|l and (|145|l are equivalent (see also Appendix [n|l . 
This yields a natural comparison |l2l | between different timescales (r for measurements and kicks, 1/K for continuous 
coupling) . 

The ratio ^'^{K)/^ is shown in Fig.|Slas a function of 2'k/K. The transition between these two regimes takes now 
place aX K = K* where K* is defined by the equation 

7^(if*)=7^(0)=7- (155) 

If K* is in the asymptotic region H153|l 

k{K*) ~ ^ = 2K'^{n). (156) 

TT 

For the exponential form factor H76|) one gets 
while for the polynomial form factor (|77|l one gets 



One observes a dependence of K* on the tail of the form factor. The strong coupling region is shown in Fig. |5| 



D. Comparison among the three control strategies 



There is a clear difference between bona fide projective measurements and the other two cases, BB kicks and 
continuous coupling. In the former case Eqs. H139|) - (|140|) yield a dependence of r* on the global features of the form 
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FIG. 9: Continuous coupling: 'y'^{K)/'y vs 2nW/K for large K. Full line: exponential form factor 17611 : dashed line: polynomial 
form factor 17711 with n = 2. Decoherence is suppressed when 'y'^{K) < 7; it is enhanced otherwise. 



K(^{u;)/g^W 

. 6 
. 4 
0.2 

q 1 d vy 



(a) 



. 6 
. 4 
0.2 



CO 



K^{u;)/g^W 



. 4 



0.2 



UJ_ 

>W 



FIG. 10: Different features of the three control methods. Form factor (polynomial, n — 2) k^{lo) (dashed line) and form factor 
modulated/multiplied by the control "response" function (full line) for: (a) pulsed measurements, Eq. 11341 . with control 
response function rsinc [{uj — Q.)t/2] (here and in the other two cases, Q. = 0.2W); (b) BB kicks, Eq. llTlTt . with control 
response function {2/tv) T,JLo 0' + k)'^ i'^ - ^ - 7(2j + 1)) + S [ui - D. + ^{2j + 1))] [see iWi and notice that the first 2- 
3 terms of the series yield an excellent approximation]; (c) continuous measurement, Eq. 115311 . with control response function 
■k[5{uj — f2 — A') + S{uj — + K)]. The gray line is a guide for the eye and interpolates {2/tt) (j + ^) k^{lo) in (b) and ^^^{lj) 
in (c). We set r = 2tx /K = 50H^~^ (a "large" value): this yields in all cases a (controlled) decay rate that is very close to that 
obtained by the Fermi Golden rule. 



factor (i.e., its integral). By contrast, Eqs. (|149|I - H152|) and H15()l) - (|158l) "pick" some particular ("on-shell") value(s). 
This important difference, due to the different features of the evolution (non-unitary in the first case, unitary in the 
latter cases), is graphically displayed in Fig. llOl and llll where the different mechanisms of control are compared. In 
Fig. ^1 r is "large" (in units of inverse bandvifidth) and the three methods yield almost no control: one essentially 
reobtains the Fermi Golden rule 7 = 2'kk^{^), although in different ways. In Fig.^] t is "small" and the effective 
lifetime is sensibly modified, although by different mechanisms. 

The three control methods are graphically compared in Figs. I12I13I The different features discussed in Figs. IIOIIII 
yield very different outputs, clearly apparent in Fig. 1131 that can be important in practical applications: decoherence 
can be more easily halted by applying BB and/or continuous coupling strategies. These two methods yield values of 
T* (or K*) that are easier to attain. However, this advantage has a price, because BB and continuous coupling yield 
a larger enhancement of decoherence for t > t* , K < K*. The two dynamical methods perform better only when 
T < r*, K > K*. This is apparent in Fig.^] We notice that a strict comparison between continuous coupling and 
the other two methods is difficult, as it would involve an analysis of numerical factors of order one in the definition 
of the relevant conversion factors between the frequency of interruptions r and the coupling K [this factor has been 
sensibly-but arbitrarily-set equal to 27r in Figs. I12fl^ see sentence after Eq. (|154ll and Appendix [C]. 
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FIG. 11: Same as in Fig. IIUI but for r = 2-k/K = ?iW~^ (a "small" value): this yields a bona fide control of the decay rate 
(in this particular situation, decoherence is enhanced in the Zeno case and suppressed in the other two cases), (a) The control 
response function rsinc^[(tj — f2)T/2] is very broad and the effective lifetime depends on the "global" features of the form factor, 
(b) For small r all the arguments of the J-functions in (I109I I tend to oo: for well-behaved form factors (like that shown in the 
figure), only the first term contributes significantly; the controlled lifetime depends on the local features of the "tail" of the 
form factor, (c) For large K the arguments of the 5-functions in Eq. 11531 tend to ±oo and the controlled lifetime depends 
again on the local features of the "tail' of the form factor. 
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FIG. 12: Comparison among the three control methods. The graphs of Figs. 0] 1^1 and [glare displayed together. BB kicks and 
continuous coupling are more effective than bona fide measurements for combatting decoherence, as the regime of "suppression" 
is reached for larger values of r and K"^ . 



VIII. SUMMARY AND CONCLUDING REMARKS 



We have analyzed and compared three control methods for combatting decoherence. The first is based on repeated 
quantum measurements (projection operators) and involves a description in terms of nonunitary processes. The 
second and third methods are both dynamical, as they can be described in terms of unitary evolutions. In all cases, 
decoherence can be halted by very rapidly/strongly driving or very frequently measuring the system state. However, 
if the frequency is not high enough or the coupling not strong enough, the controls may accelerate the decoherence 
process and deteriorate the performance of the quantum state manipulation. The acceleration of decoherence is 
analogous to the inverse Zeno effect, namely the acceleration of the decay of an unstable state due to frequent 
measurements [43Ll43|. 

As a general rule, when one endeavors to control decoherence by suitably tailoring the coupling of the system 
of interest to another system (such as an external field, or a measuring apparatus), one should carefully look at 
the relevant timescales, as it is not true that repeated measurements/interruptions always lead to a suppression of 
decoherence. 
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FIG. 13: Comparison among the three control methods: small times/strong coupling regions. The graphs of Figs. Eland H 
are displayed together. 



It is convenient to summarize the main results obtained in this article in the particular case of a two-level system 
(qubit) with energy difference fl. If the frequency of measurements or BB kicks, or the strength K of the coupling 
tend to oo, the two-dimensional (Zeno) subspace defining the qubit becomes isolated and decoherence is completely 
suppressed. However, if and K are large, but not extremely large, the transition (decay) rates between the 
qubit subspace and the remaining sector of the Hilbert space display a complicated dependence on r^^ and K, and 
decoherence can be suppressed or enhanced, depending on the situation. 

At low temperature (3~^ ^ fl ^ W , where W is the bandwidth of the form factor of the interaction, the decay 
rates read, from (P^) . ((n5|) . fl^ . fTIS)) and fTC^ 

^ z 

< jHr)-^n{^), r^O, (159) 
^ -/%K) - 7rn{K), K ^oo , 

where Z, k and c denote (Zeno) measurements, (BB) kicks and continuous coupling, respectively, n is the form factor 
and l/r^ — / dujK,{LL>) the Zeno time (more accurate definitions were given in the preceding sections). As we have 
shown, there is a characteristic transition time r* [coupling K*], such that one obtains: 

decoherence suppression : 7(t) < 7 [■y{K) < 7], for t < t* [K > K*], 
decoherence enhancement : 7(t) > 7 [^{K) > 7], for t > t* [K < K*]. (160) 

Therefore, in order to obtain a suppression of decoherence, the interruptions/coupling must be very frequent/strong. 
Notice, in this context, that both t* and 2tt/K* are not simply related to the inverse bandwidth 2'kW~^: they can 
be in general (much) shorter. For instance, in the Ohmic polynomial case (|77ll . one easily gets from (|78f) and H159|l 

' T*z ~ 2^W-^ {2{n ~ l)al§) « 2^W-^ , 

, r*c,2^W-^^[^^)^' <^2nW-\ (161) 

where a„ = (\/7r/2)r(n — 3/2)/r(n — 1) < 7r/2 is a coefficient of order 1 and n characterizes the polynomial fall off 
of the form factor (|77|) . The above times/coupling may be (very) difficult to achieve in practice. In fact, we see here 
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that the relevant timescale is not simply the inverse bandwidth 2ttW~^, but can be much shorter if <C W, as is 
typically the case. These conclusions, summarized here for the simple case of a qubit, are valid in general, when one 
aims at protecting from decoherence an Ai"-dimensional Hilbert subspace. 

An important example that we have not explicitly analyzed in this article is the case of 1// noise, and its sup- 
pression by means of techniques like those discussed here. There has recently been a surge of interest in this issue in 
quantum information processing devices, where such noise is often attributable to (but certainly not limited to) charge 
fluctuations in electrodes providing control voltages [67l l6^ . The need for such electrodes is widespread i n q uantum 
computer prop osals, e.g., trapped ions (where observed 1// noise was reported in 69]), quantum dots [73, doped 
silicon [lEIlil, electrons on hehum [t^, and superconducting qubits ^3- latter case, in a recent experiment 

involving a charge qubit in a small superconducting electrode (Cooper-pair box), a spin-echo- type version of BB was 
successfully used to suppress low- frequency energy-level fluctuations (causing dephasing) due to 1// c harp e noise . 
Several recent papers have dealt with suppression of this particular kind of noise via BB decoupling 41, 64 1 |75l l76l| . 
The "bottom-up" approach models 1// noise as arising from a collection of bi-stable fluctuators ,41j, ,67^ .751 l7a | . The 
alternative is to treat 1// noise as contributing a particular form factor 0) El| ■ We will pursue these ideas as a 
future topic of investigation, but we expect that the main results obtained in the present paper be applicable to this 
case as well. 

The results obtained in this paper are of general validity and bring to light the different features of the control 
procedures as well as the crucial role played by the form factor of the interaction. We do not expect any drastic 
change for different decoherence mechanisms and/or different physical systems. The only somewhat delicate issue, 
in our opinion, is to understand whether the system investigated can be consistently described by means of a set of 
discrete levels. 
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APPENDIX A 



In this Appendix, the assumption of the factorized form of the initial density operator, as in jSJ, which is usually 
taken for granted, is shown to be justified in the weak-coupling (scaling) limit. We only outline the main derivation. 
Further details will be reported elsewhere [i^. 

Consider the initial-value problem 

d 

—p = CtotP = {Co + \Csb)p = {Cs + Cb + XCsb)p, p(0) = po, (Al) 

where the dependence on the coupling constant A of the interaction Liouvillian Csb is made explicit. Notice that 
the initial density operator can be of any form and is not assumed here to be factorized like in Q. The projection 
operators V and Q, defined in (|41|l . and the above Liouvillians Cs, Cb and Csb satisfy the same conditions (|42|) and 
H43|l. The projected density operators Vp and Qp satisfy 

d d 

—Vp = CoVp + \VCsbQp, ^^Qp^{Co + \QCsbQ)Qp + ^QCsbVp, (A2) 
ot at 



respectively. Following the same procedure as in Sec. lIIII we arrive at the following exact equation for the 7^-projected 
operator in the interaction picture 



d 



-Cot 



dt 

( dt'e-^°'VCsBe'^^°+^^^'"^^''CsBVpit~t'). (A3) 
Jo 
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Notice that the first term on the rhs represents the contribution arising from a possible initial correlation between 
the system and reservoir. We now show that this term dies out in the weak-coupling (i.e., scaling) limit A — > with 
fixed T = X^t. For this purpose, define 



pi(T;A)^e-^«-/^Vp(r/A2), 



(A4) 



that satisfies 



pi(r;A) = ie-^''-AV/:sBe(^''+^e£.BQ)r/A^g^^ 
A 

+ / dt'e-^«^/^V/:si3e(^°+^2^^-2)*'/:5Be^°(^/^'-*')7'pi(T - \H'; A). 

^0 

The first term vanishes in the A ^ limit 1461 , since 



/•OO -1 /"CXD 

/ dT-e^^/^'Y{T) = A / dTe^^Y{\^T) ^ 0, 
Jo Jo 



as A ^ 0, 



(A5) 



(A6) 



for any superoperator such that the integral Jq°° dre'^'^ exists. This means that the contribution originating from the 
initial correlation between the system and reservoir disappears in the scaling limit and therefore we are allowed to 
start from an initial density matrix in the factorized form 
Finally, the dynamics of pi(t; 0) is governed by 



pi(r; 0) = V / dt'VCsB{0)QCsBht' 
Jo 



)Q^Pi(t;0) 



(A7) 



with the factorized initial condition 101, where the Q^^ are the eigenprojections of the Liouvillian Cs defined in H52|l . 

From a physical point of view, the factorization ansatz described in this appendix simply means that the "initial" 
correlations between the system and its environment are "forgotten" on a time scale of order A^. We also note that 
several authors have addressed the question of the modifications that arise when it is not permissible to assume 
initially separable system-environment, e.g., (t^ . 



APPENDIX B 



We derive Eq. (|57|l . The first equality reads 



p^Ctatrp 



t/r 



PVz{t)P 



Let us write Vz{t) — F(r, r), where 



Vz{t) = e^o^Texp ( j ds e-^<>''gz{T)e 



CnS 



V{t,u) = e^''*Texp (^J^ ds e-^"'gz{u)e^<'' 



By deriving with respect to i, we get 

dtV{t,u) = [Co+gziu)]V{t,u), 

so that 

V{t,u) ^ eiip{[Co + GziuM, 
where we used ^(0,-^) = 1. As a consequence, Vz{t) ~ exp{[£o + gz{T)]T} and 

PgAotrp] ^ \p g^p| ^ {t)]t}p] = P exp{ [Co + Qz (t)]0, 



(Bl) 

(B2) 

(B3) 
(B4) 

(B5) 



because \P,lio\ = [P.Gz] = 0. This is Eq. 
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Let us now solve Eq. 

f (ii e~^^*/:z(T)e^^* = / dtPCi{t)P = f dt f dsPJCi{t,s)P. 
Jo Jo Jo Jo 

By using (gHjl and (O 

ICi{t, s)a = tYB {CsB{t)CsB{s) <7<E>pb}, CsB{t) = e-^^'+^^'^' CsBe^^'+^^'^\ 

we get 

^ e-^s*£z(T)e^^*a = ^ dt ds Iyb [PCsB{t)CsB{s) Pa ® pe] 

= dte'^^'J dstrB{PCsBCsB{s)Pe^''a(g)pB} ■ 
Let us rewrite the previous equation in terms of the eigenprojections of Cs defined by H52|l : 

J2f dte'^^-^'^'QMr)Qu'a = Y, j dte'("-'"')* / ds tvB [Q^PCsBCsB{s)PQ^'<y ® pb] 
Performing the first integral, we get 



(B6) 
(B7) 



(B8) 



(B9) 



(BIO) 



(Bll) 



ds tr_B |Q(^^>CsB'C5B(s)F(3i^'Cr (g) psj , 5(2:) 

Since g(0) = 1, the diagonal terms yield 

QuCz{t)QuO = ^ j dt j ds tiB \^QujPCsbCsb{s)PQujO- pB^ ■ 

The off-diagonal terms do not contribute to the master equation, as explained at the end of Sec. IIIII Eqs. (|50|I - H55() . 
By using the property and noting that [P, Q^] = by fT^ . we get 

£z{t)<j = -J2- f dt f dstrB{p[[Q^HsB) ,[[q^u,Hsb{s)) ,P<j<»pb]]}, 

whence, by using ()56|l . 



(B12) 



E 



dt I ds trB \ P 

dt I ds tiB |-P 
+P 



H 



(m) 
SB ' 



H^f/'\s),Pa<g>pB 



} 



Xm® A]^{s),Pa® pB 



(B13) 



where, like in Eq. (|67|l . in the second equality we neglected terms containing two annihilation or creation operators. 
From (|BT3|) we get Eq. ((HU with 



By noticing that 



we finally get 



7^(r) = - Re dt I ds [{A„MAl{s)) + (aL™(0)A_,„(s))) 

du; K^(u;)e*(---'")^ 



7™(r) = 
which is Eq. (|^ of the text. 



(A^Alis)) + (^Al^A_ra{s) 

2 f°° g I- cos(w - w,„)r 
duj 

, {U}-U!rnY 



dw nti^) 



sm {-^ 



(B14) 



(B15) 



(B16) 
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APPENDIX C 

We derive Eqs. ((TIIT|l and ifTU^ . We start from Eqs. (^3 and ^ 

f rfs e-^-^^k(r)e^-^ = / ds Csb{s) , 
Jo JO 



ds e-^-"gk(r)e 



ds / dsi 



Csb{s)Csb{si) - e-'^-'T\,{T)e 



(CI) 
(C2) 



where Ct = ^ + Cq, and by taking the trace over the bath we get 



dt e'^s*'Ck{T)e^'s*'a = / dt f ds trs 



^ dt ds tYB { [e-^'^'CsBCsBis) - e-<''*J^k(r)e-^-^J^k(r)e^-^e^'s*] a ® ps} , (C3) 



with 



T 



(C4) 



Equation (|C3|I is similar to (|B8p and, by projecting onto the eigenprojections Puj{t) of Cs{t) and taking only the 
diagonal terms, one obtains Eq. H101|l . However, in order to compute the decay rates 7m(T) one can give an alternative, 
more physical derivation by elaborating on the technique of Ref. First notice that the BB dynamics H93|l is 

generated by the time-dependent Hamiltonian 

H (t/T) = Htot + Sp it/T) , (5p (a;) = ^ 5(x - n) . (C5) 

nez 

In the enlarged Hilbert space TL ® L?{T) we can consider the (time- independent) Floquet Hamiltonian 



where 
We get 

whence S "H, 



T T 

0e [-1/2,1/2), pe = -ide, [0,Pe]=i- 

= -t[0,HF,oq] = l/T, 0{t)=t/T, 



(C6) 

(C7) 
(C8) 



A{t) = -i [A{t), Fpioq] = -i [A{t), H (t/T)] , (C9) 

so that every observable in 7i evolves according to the original Hamiltonian (|C5|) . The eigenvalue equation for pe 
reads 



pe\m) = 27rm|m), 



(CIO) 



The Hamiltonian ljC6|l in 7i ® L'^(T) represents a control by a strong continuous coupling, analogous to that discussed 
in Sec. I VII if one identifies K ~ 1/t and He = pe- Therefore, from Eq. p22|l and (|C10|I we obtain 



,(t) = + + o(r) = ^ + f7W + 0(r) , 



and from Eq. H119() we get 



7L.(r) = 2n4 (c.„,„(r)) = 2n4 + + 0(r)) 



(Cll) 



(C12) 



which is Eq. (|103|l of the text. 
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